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ABSTRACT 

A small fraction of the halo field is made up of stars that share the light element (Z < 13) anomalies characteristic of second generation 
globular cluster (GC) stars. The ejected stars shed light on the formation of the Galactic halo by tracing the dynamical history of the 
clusters, which are believed to have once been more massive. Some of these ejected stars are expected to show strong A1 enhancement 
at the expense of shortage of Mg, but until now no such star has been found. We search for outliers in the Mg and A1 abundances 
of the few hundreds of halo field stars observed in the first eighteen months of the Gaia-ESO public spectroscopic survey. One 
halo star at the base of the red giant branch, here referred to as 22593757-4648029 is found to have [Mg/Fe] = -0.36 ± 0.04 and 
[Al/Fe] = 0.99 ± 0.08, which is compatible with the most extreme ratios detected in GCs so far. We compare the orbit of 22593757- 
4648029 to GCs of similar metallicity and find it unlikely that this star has been tidally stripped with low ejection velocity from any 
of the clusters. However, both chemical and kinematic arguments render it plausible that the star has been ejected at high velocity 
from the anomalous GC m Centauri within the last few billion years. We cannot rule out other progenitor GCs, because some may 
have disrupted fully, and the abundance and orbital data are inadequate for many of those that are still intact. 

Key words. Stars: abundances - Stars: Population II - Techniques: spectroscopic - Globular clusters: general - Galaxy: halo - 
Galaxy: stellar content 


1. Introduction 


The majority of stars in the Galactic halo reside in the field, 
leaving only 2% of the total stellar mass in globular clusters 
(GCs, see e.g. Vesperini et al. 2010 and references therein). 


* Based on data acquired by the Gaia-ESO Survey, programme ID 
188.B-3002. Observations were made with ESO Telescopes at the La 
Silla Paranal Observatory. 


How large a fraction of today’s field stars escaped from a clus¬ 
ter and the original sizes of these clusters is a matter of much 
debate. With few exceptions, field and cluster stars share similar 
chemical compositions for alpha, Fe-peak and neutron-capture 
elements ( |Gratton et al.| 12004) [Pritzl et al.||2005| l. However, the 
GCs display unique abundance patterns for lighter elements that 
are not typically seen in the field, in particular for C, N, O , Na, 
Mg, and Al. The existence of such anomalies have been known 
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for decades ( [Cottrell & Da Costa| |1981[ |Kraft et aT] |1997|) and 


have been mapped in detail in large GC samples (e.g. Carretta 
et al. 2009 |l. We now know that multiple populations of stars re 


side within each cluster, with some stars having been polluted by 
the elements produced during H-burning at high temperatures by 
more massi ve stars (| Denisenkov & Denisenkova|[T989||Ventura| 
|et al.| |2001t [Decressin et al.f|2007 1. Consequently, a large frac¬ 
tion of present-day GC stars are enhanced in N and Na, at the 
cost of depleted levels of C and O. The second generation stars 
are believed to also be enhanced in He, the main product of H 
burning. 

The correlating and anti-correlating patterns of the lighter 
elements characteristic of GCs are seen neither in disk open 
clusters ( |de Silva et al. 2009[ Bragaglia et al. 20141 nor in 
the Galactic bulge ( [Lecureur et al. 2007 [ Bensby et al. |2013| l. 
Therefore, by finding out how many field stars have light element 
anomalies and thus likely originated in GCs, we can study the 
GC-field link. The largest systematic studies undertaken so far 
are those by Martell & Grebel P010| l and Martell et al.| ( |20lT] l, 
who focussed on C and N anomalies as traced by molecular 
bands in SDSS/Segue spectra. The latter study reports that 3% 
of field stars display the chemical characteristics of GC stars. 
This small number agrees well with higher-resolution studies 
with fewer number statistics, e.g. Carretta et al. (2010; 1.4%) 
and [Ramifez et ak (2012; 3 + 2%). The Ramirez et al. study 
could confirm both elevated Na and depleted O-abundances in 
two halo field dwarfs, a clear indication of a GC connection. 
As discussed, e.g., by |Gratton et al. ( 2012| l, models of halo for¬ 
mation must account for much higher initial masses of GCs to 
explain this fraction of second generation GC stars in the field 
(however, see Prantzos & Charbonnel 2006 and Bastian et al. 
2013). 

The O-Na anti-correlation was found in all GCs where it has 
been investigated with enough statistics, indicating that the pro¬ 
genitor objects were massive enough to activate the NeNa-cycle. 
However, the anti-correlating abundance pattern of Mg and Al 
that is expected from the operation of the MgAl-chain manifests 
itself less uniformly. Several clusters display a small Mg-spread 
and a significant Al-spread ( [Carretta et al.| 2009 j l while very few 
clusters, e.g. cu Cen and NGC 2808, show a pronounced Mg-Al 
anti-correlation (see Fig.|^. The likelihood of finding escaped 
stars in the field with clear GC characteristics for both Mg and 
Al is thus small and it requires a large spectroscopic survey to 
detect them in significant numbers. Here we report on the first 
discovery of a halo field star with both strongly depleted Mg and 
strongly elevated Al abundance. 


2. Observations and analysis 

The Gaia-ESO public spectroscopic survey ( [Gilmore et al. 
2012| l is a five-year survey with FLAMES on the VET, which has 
been observing since December 2011. Two multi-object spectro¬ 
graphs operate simultaneously in high resolution (UVES) and 
medium resolution (GIRAFFE/MEDUSA), with fainter targets 
such as distant halo stars, allocated to the GIRAFFE fibres. Two 
settings are used for the halo stars; HRIO, which covers 535- 
565 nm at a spectral resolving power of R - A/6A - 19 800, 
and HR21, which covers 845-900 nm at R - 16 200. The com¬ 
bined analysis of data in the two GIRAFFE settings allows for a 
determination of stellar parameters, Al, alpha-elements, and Fe- 
peak elements. For 22593757-4648029, we achieve an average 
S/N per pixel of 25 in HRIO and 55 in HR21 after 2x1500 s of 
exposure in each setting (see sample spectra in Fig.[^. The ob¬ 
servational data are summarised in Table 1. A full account of the 


Table 1. Data for 22593757-4648029. 


Name 

Value 

Name 

Value® 

RA 

22 59 37.57 

Teff 

5261 ±36± lOOK 

DEC 

-46 48 02.9 

log^ 

2.84 ± 0.54 ± 0.25 

k 

15.649 ±0.037® 

[Fe/H] 

-1.49 ±0.05 ±0.20 

K 

13.928 ±0.061® 

[Mg/Fe] 

-0.36 ±0.04 ±0.10 

J-K 

0.494® 

[Al/Fe] 

0.99 ±0.08 ±0.10 

TcSJ-K 

5326 ± 140 K® 

[Si/Fe] 

0.34 ±0.13 ±0.10 

Eb-v 

0.010® 

[Ca/Fe] 

0.41 ±0.13 

krad 

44.97 ±0.16kms-'® 

[Ti/Fe] 

0.46 ± 0.22 

Pra 

2.8 ± 1.8mas/yr® 

[Y/Fe] 

< 0.11 

Pdec 

-10.40 ± 1.8 mas/yr® 

Distance 

6.7+1® kpc 


Gaia-ESO Survey 

(d) 

Casasrande et al. 1 2010 

1 

f*” APASS 

(e) 

Schlegel et al. 



2MASS, Skrutskie et al. (2006 


target selection, data reduction and processing will be given in a 
dedicated GIRAFFE data release paper. 

Several independent analysis nodes, using independent 
methods, participate in the determination of stellar parameters 
and chemical abundances in the Gaia-ESO Survey. The re¬ 
sults of individual nodes are thoroughly examined, evaluated, 
and homogenised to arrive at the finally recommended values 
for each star. We adopt the homogenised, recommended re¬ 
sults for 22593757-4648029, based on spectra from the inter¬ 
nal data release 2 (GESviDR2Final). The object is classified 
as a metal-poor star at the base of the red giant branch with 
Teff = 5250 + 34 + 100K, logg = 2.84 + 0.54 + 0.25 and 
[Fe/H] = -1.49 + 0.05 + 0.20 (random and systematic error) 
. The details of all the analysis methods and the homogenisa¬ 
tion procedure will be described in Recio-Blanco et al. (in prep). 
The distance to the star was determined based on Bayesian fits to 
evolutionary tracks using our derived spectroscopic parameters 
following the methodology described in[Serenelli et al.[(|2013|l. 


3. Discussion 


The peculiar chemical composition of 22593757-4648029 is 
reported in Table While three of the alpha-elements. Si, 
Ca, and Ti, are compatible with the standard 0.4 dex [a/Fe]- 
enhancement characteristic of metal-poor stars. Mg falls approx¬ 
imately 0.8 dex below the expected value. In contrast, Al is in¬ 
stead enhanced to 1 dex above solar. As discussed below, this is 
compatible with the most extreme populations found in GCs and 
we argue that the star is a GC escapee. It was discovered in a 
sub-sample of ~ 7300 FGK stars (half of all GIRAFFE targets) 
with detectable Mg and Al abundances in the disk and halo field. 
Based on a simple two-component fit to the metallicity distribu¬ 
tion function of this sub-sample (peaks at [Fe/H] = -1.6 + 0.5 
and [Fe/H] = 0.24 + 0.36) , we expect an approximate halo 
fraction of 5-10%, i.e. a few hundred stars. Finding one such 
chemically unusual star among them is not inconsistent with es¬ 
timates of GC escapees in the halo of 3% (Martell et al. 2011 1 , 
noting that there are also two other marginal candidates with 
[Mg/Fe] ~ 0, [Al/Fe] ~ 1, and [Fe/H] = -0.6/ - 1.3. We omit 
these additional stars from the discussion and plots because of 
the significantly lower S/N of their spectra. 

While we consider a GC origin of 22593757-4648029 the 
most likely explanation for its non-typical abundance patterns, 
one may speculate on alternative reasons, e.g., mass transfer 
from a binary companion in the field. The so-called CEMP- 
s stars (e.g. [Lucatello et al.[ [2005|l are a class of metal-poor 
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Wavelength [A] 




8771 


8772 8773 8774 8775 

Wavelength [A] 


8776 


Fig. 1. Observed (dotted line) and best-fit synthetic (solid line) spectra of 22593757-4648029. Shaded regions indicate +0.3 dex in 
abundance. 


stars believed to have been polluted with gas transferred from 
an AGB companion. These are characterised by enhancement in 
carbon and slow neutron-capture elements. The spectral range 
does not allow a stringent constraint on [C/Fe], but the upper 
limit [Y/Fe] < 0.11 excludes strong s-process enhancement. 
Further, CEMP stars are not characterised by low Mg abun¬ 


dances; no star in the sample presented by Allen et al. (2012 


including literature studies) has sub-solar [Mg/Fe]. According to 
Ventura et al. (20111, AGB stars with ~ 6Mo produce the most 


extreme Mg-Al-Si nucleosynthesis. If such a star was once a bi¬ 
nary companion to our presumably old and low-mass field star, 
the system would have an unusual mass ratio. 

Low [Mg/Fe]-ratios are commonly found in dwarf 
spheroidal galaxies at this metallicity, but are then accom¬ 
panied by similarly low ratios of other a elements, like Ca, Si, 
and Ti, with respect to Fe (see e.g. |Koch & McWillial^ [2008 [ 
|Kirby et al.[ |2009 1 . No [Al/Fe] enhancement has been found 
(nor is it expected) in these systems. 

In Fig.|^ we compare the Mg and Al abundances of the 
field star to six different GCs. M22 and NGC36^ like most 
GCs, display no striking Mg-Al anti-correlation and the lowest 
[Mg/Fe] values of these clusters are far from that of 22593757- 
4648029. NGC6752, NGC2808 and co Cen all show strong ev¬ 
idence of multiplicity and internal He variations, having at least 
triple main sequences and extended horizontal branches (e.g. 

Milone et al. 2013[ King et ak] 2012| l. As 


Piotto et al.[ 12007 




seen in Fig.l^ the [Al/Fe] ratios in the most highly polluted stars 
in these clusters are as high as that of 22593757-4648029, while 
the [Mg/Fe] ratios of stars in NGC 2808 and a> Cen are almost 
as low. However, NGC 2808 is too metal-rich to well match the 
metallicity of 22593757-4648029, making oj Cen the more likely 
parent GC. Unfortunately, there are surprisingly few Mg abun¬ 
dance measurements published for this otherwise well-studied 
GC. The fact that [Mg/Fe] ratios as low as -0.36 can be found in 
Galactic clusters is evidenced by the peculiar object NGC2419, 
where the record holder has [Mg/Fe] a; -1 ( jMucciarelli et al.] 
2012 Cohen & Kirby[ 2012| l. This cluster, however, is instead 


http://WWW.astro.yale.edu/dana/gc.html 
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too metal-poor to be a plausible parent. 

It is easily realised that chemistry alone in not conclusive, 
because as many as 50 GCs have [Fe/H] within +0.20 dex of 
the field star ( |Harris[[l996[ 2010 edition) and most of them lack 
Mg and Al data. However, further insight can be obtained from 
the known kinematics of the star (see Table[T]l and about half of 
the candidate GCs with matching metallicity (D. CasettQl. We 
proceed under the assumption that the progenitor cluster is still 


Fig. 2. Mg and Al abundances of normal field stars and five 
GCs compared to 22593757-4648029 (red square with error 
bars). The average metallicity is listed after the cluster name. 
Abundance data represented with plus signs are taken from 
[Worley & Cottrell|(2010; NGC 362),|Marino et al.|(|2009l 2011; 
----- - - . ■ ^, 2005 ; 

(2009; 


M22), Norris & Da Costa (1995; at Cen), Yong et al. 
NGC 6752), jCarrettal (2006; NGC 2808), jCarretta et alT 


* Mg and Al abundance data for NGC 362 were determined using 
spectra and stellar parameters obtained byjWorley & Cottrell j2010| . 


NGC 2808), and [Cohen & Kirby| (2012; NGC 2419). GES rec¬ 
ommended data are marked with bullets. Arrows indicate upper 
limits. 


intact, while evidence has also been found in the inner halo for 
disrupted GCs ( [Bernard et al. 2014| l. 

We assess the likelihood that the star was previously inside a 
known GC by integrating the current orbit of the star and the GC 
5 Gyr back in time in the Milky Way potential and look for close 
encounters. Eollowing the work of e.g. jJohnston et al.| ( 1999| l and 
[Price-Whelan & Johnston[(|2013[), the energy of the star during a 
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close encounter is computed as £ = where 6V 

is the relative velocity of the star to the cluster and 6r is the dis¬ 
tance from the cluster to the star. Encounters within rcc and with 
£ < 0 are consistent with the star having been tidally stripped 
with a relative velocity smaller than the escape velocity, while 
encounters within rcc and £ > 0 are consistent with the star 
having been ejected with a velocity of approximately V2 £. We 
proceed to Monte-Carlo sample the proper motions, radial ve¬ 
locities, and distances of the star and of the clusters according 
to their error bars. Then, for each Monte-Carlo sample, we de¬ 
termine whether encounters occurred within the tidal radius and 
record the star’s energy at the time. The fraction of Monte-Carlo 
samples with close encounters and £ < 5 therefore gives 

us a probability that a star was ejected from a given cluster with 
velocity V < Vgection i-e. £(ejected,V < Vgectionl GC). The plot 
of this cumulative probability is shown in Figure]^ 


Evidently, given the current observational errors, it seems 
unlikely that the star has been tidally stripped from any clus¬ 
ter i.e. with small ejection velocity. But if we allow it to have 
been ejected at high velocity (> lOOkm/s), which is possible 
in some scenarios involving black holes and binary system in¬ 
teractions (Gvaramadze et al. 2009 [ Liitzgendorf et al. 2012| l, 
a handful of clusters could possibly have been associated with 
22593757-4648029 in the past. Normalising the total probability 
to one, the top three posterior probabilities of ejections assum¬ 
ing ^ejection <200 km/s is NGC5139 (w Cen) 44%, NGC6656 
(M22) 12%, NGC 362 12% while for some clusters this proba¬ 
bility is negligible, e.g. NGC 5272, NGC 5904, and NGC 7006. 
More details of the orbit calculations are listed in Appendix A. 

Interestingly, we note that a> Cen, which appeared most plau¬ 
sible already from the limited chemical information at hand is 
also favoured based on pure dynamical arguments. A more de¬ 
tailed chemical abundance analysis, including more elements, 
would help to verify this hypothesis. In particular, Mn and Cu 
abundances may help because a fraction of ai Cen stars are 
under-abundant in these elements compared to the held and 
other GCs ( |Cunha et al.| |2010[ [Romano & Matteucci| [2007 j 
and references therein). In a few years time, the Gaia mission 
will have greatly improved the accuracy and statistics of dis¬ 
tance and proper-motion information for stars and star clusters. 
Combining this information with radial velocities and chemistry 
from ground-based telescopes, we will be able to connect stars 
like 22593757-4648029 to their parent GCs with greater conh- 
dence. 
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Appendix A: Orbit caicuiations 

The Milky Way potential that we used in the simulations is the 
analytical 3-component bulge-, disk-, and halo-potential used to 
model the orbits of satellites in [Johnston et ^ ( |1995| l, |Dinescu| 
et al. (19991, and |Koposov et akf i 2010) 1. The potential of each 
cluster was approximated by a Plummer sphere using currently 
measured half-light radii and masses. The calculation of the tidal 
radii was done assuming constant rotation curve (e.g Milky Way 
enclosed mass is proportional to radius). 

Since the true Galactic potential is not accurately known and 
the results of our simulations may depend significantly on the 
assumed parameters, we also run the simulations with masses 
of the disk and bulge perturbed by 10% Gaussian variations and 
see how it affects the probabilities of the star being associated 
with one of the clusters. As an additional test, we also run our 
backwards orbit integration to 2 Gyr instead of 5 Gyr to see how 
it affects our results. While the individual probabilities for the 
clusters do change, the overall picture does not, e.g. u Cen is 
always ranked within the top three. 

We further note some important limitations of our calcula¬ 
tion: we ignore the dynamical friction of the GCs, which is ex¬ 
pected to bring the GCs closer to the Galactic centre. We also ig¬ 
nore the fact that all the GCs were likely to be bigger in the past 
and have tidally lost some of that mass, which reduces the abso¬ 
lute probabilities we find. Finally, we ignore secular changes in 
the Milky Way potential over time, which are expected although 
the Galaxy had a quiet recent accretion history. 
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